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I .  INTRODUCTION 


In  the  Fourier  integrals  used  to  calculate  stresses  in  a  hollow 
cylinder,  the  variable  of  integration  ranges  from  zero  to  infinity  along 
the  real  axis;  consequently  the  Bessel  functions  in  the  integrand  must 
be  computed  for  very  large  and  very  small  values  of  the  argument. 

Calculation  of  stresses  in  a  thin-walled  cylinder  by  residue  theory 
requires  Bessel  functions  of  complex  argument  and  very  large  modulus. 

Use  of  the  Bessel  function  previously  described1 2  to  calculate  the  integrand 
for  these  extreme  ranges  leads  to  exponential  overrun  and  cancellation 
of  significant  digits.  Axial  stresses  due  to  discontinuous  internal 
shear  loading  are  calculated  from  a  Fourier  integral  with  an  apparent 
singularity  at  the  lower  limit  of  integration.  Evaluation  of  the  integral 
by  quadratures  leads  to  the  numerical  difference  of  nearly  equal  numbers 
for  points  near  the  discontinuity  of  loading. 

Exponential  overrun  is  especially  critical  in  calculating  stresses 
for  thin-walled  cylinders  according  to  the  theory  of  elasticity. 

Asymptotic  series  have  been  used  to  calculate  stresses  for  cylinders  with 
wall  ratios  in  the  range  1.01  <_  B  <_  1.25.  The  lower  part  of  this  range 

is  not  accessible  when  conventional  tables  or  subroutines  are  used  to 
calculate  the  integrand  of  the  Fourier  integral. 

This  report  describes  three  algorithms  which  have  been  used  to  over¬ 
come  the  numerical  difficulties  described  above.  The  corresponding  computer 
programs  are  listed  in  the  thick-walled  cylinder  code^  and  will  not  be 
repeated. 

II.  SOURCES  OF  NUMERICAL  DIFFICULTIES 

The  characteristic  function  for  a  hollow  cylinder  with  stress-free 
cylindrical  surfaces  illustrates  the  numerical  difficulties  in  evaluating 
determinants  arising  in  thick-walled  cylinder  theory. 


Vp) 

k0Cp) 

a^Cp) 

a1K1 (p) 

ijCp) 

“Kj (P) 

P^Cp) 

-PK0 (p) 

K0(q) 

ViCq) 

BjKj (q) 

ij(q) 

-Kj(q) 

qi0(q) 

-qK0(q) 

where 

1K.L.  Zimmerman,  A.S.  Elder ,  A.K.  Depue,  "User's  Manual  for  the  BRL 
Subroutine  to  Calculate  Bessel  Functions  of  Integral  Order  and  Complex 
Argument",  Ballistic  Research  Laboratory  Report  No.  ARBRLTR-02068 
May  1978.  (AD  A056369) 

2 

A.S.  Elder  and  K.L.  Zimmerman,  "Stresses  in  a  Gun  Tube  Produced  by 
Internal  Pressure  and  Shear",  BRL  Memorandum  Report  No.  2495,  June 
1975.  AD  M012765. 
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and 


P  =  sa  , 

q  =  sb  , 

a  =  p  +  (2-2v)/p  , 

=  q  +  (2-2v)/q  , 


v  =  Poisson !s  ratio  , 
a  =  inner  radius  of  cylinder  , 
b  =  outer  radius  of  cylinder. 

If  we  evaluate  A(s)  by  approximating  the  Bessel  functions  with  Hankel 
asymptotic  series,  we  find  products  of  the  type 


sa  -sb 
e  e 


and  e 


-sa 


e 


sb 


occurring,  leading  to  exponential  overrun  when  |sa|  and  | sb j  are  very 
large.  We  also  find,  after  considerable  algebra,  that  the  first  and 
second  terms  of  the  Hankel  asymptotic  series  are  lost  by  subtraction, 
leading  to  a  loss  of  significant  figures. 

The  characteristic  function  A(s)  has  a  pole  at  the  origin,  but  is 
analytic  elsewhere  in  the  finite  part  of  the  complex  s  plane.  However, 
the  Bessel  functions  of  the  second  kind  in  Eq.  (1)  have  logarithmic 
singularities  at  the  origin.  The  logarithmic  singularities  will  cause  a 
loss  of  significant  figures  when  |s|  is  small.  The  algorithm  for  calculating 
cross  products  of  the  modified  Bessel  functions  for  small  |s|  eliminates 
this  difficulty. 

III.  ASCENDING  SERIES  FOR  CROSS  PRODUCTS  OF  MODIFIED  BESSEL  FUNCTIONS 

Cross  products  For  the  modified  Bessel  functions  are  defined  by  the 
following  equations^ 


o,o(s:i  * 

I0(sa)K0(sb)-K0(sa)I0(sb) 

I0(sa)K1(sb)+K0(sa)I1(sb) 

l,0(s)  = 

IjCsa)  K0(sb)+K1(sa)I0(sb) 

1,1 

I1(sa)K1(sb)-K1(sa)I1(sb)  , 

3 

S .  Timoshenko  and  J.N.  Goodier,  Theory  of  Elasticity ,  Second  Edition 3 
McGraw-Hill  Book  Company 3  New  Iork3  1951 . 
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where  s  is  the  variable  of  integration  in  the  Fourier  integrals.  We  now 

show  that  the  L  n(s)  functions  are  analytic  in  s  except  for  a  simple 

1  ,U  4 

pole  at  the  origin.  We  have 


(sa)2  (sa)4  ^  (sa)6 


o'-  J 

22 (1 ! )  2  24 (2  ! ) 2  26 (3 ! ) 

I0(sb)  = 

(Sb)2  (Sb)4  (sb) 

1  2^24  6 

2  (1!)  2(2!)  2(3!) 

IjCsa)  = 

(sa)  (sa)3  (sa)5 

2  „3,  ,  „  5 

2  1!  2!  2  2! 3! 

I^sb)  = 

(sb)  (Sb)3  (Sb)5 

2  ,5,.,.  • 

+  . 


2  1 ! 2 !  2  21 3! 


(3) 


for  modified  Bessel  functions  of  the  first  kind.  In  the  series  for 
modified  Bessel  functions  of  the  second  kind,  we  separate  the  logarithmic 
terms  to  facilitate  the  algebraic  manipulation.  For  the  outside  radius 
we  have^ 


KQ(sb)  = 

G0(sb)  -  M0(sb)  -  S0(b,s)  , 

(4) 

where 

G0(sb)  = 

2  777  (*£ )  (:  *  r  t+  •••  *  ?)  - 

r=l  1  J 

(5) 

MQ (sb)  = 

I0«sb’  10h  (?)  " 

(6) 

S0(b,s)  * 

!  I0(sb)  [y  +  loge  s]  , 

(7) 

4 

V. K.  Prokopov,  "Equilibrium  of  an  Elastic  Axisymmetrically  Loaded  Thick- 
Walled  Cylinder ",  Priladnaya  matematika  i  mehanika ,  Vol.  XIII,  1949- , 
pages  135-144.  Institute  of  Mechanics  of  the  Academy  of  Sciences ,  USSR. 
FTIO  Translation  No.  1-2589,  Aberdeen  Proving  Ground ,  Maryland, 
Translation  dated  22  August  1967. 

** British  Association  for  the  Advancement  of  Science,  Bessel  Functions, 
Part  I,  Function  of  Orders  Zero  and  Units,  Mathematical  Tables,  Volume 

VI,  University  Press,  Cambridge,  1937. 
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where  y  is  Euler !s  constant. 


Similarily  we  have 


(sb)  =  -  G^sb)  +  M1(sb)  +  S^(b,s)  , 


where 


Gj (sb) 


(sb)  =  Ij(sb)  loge  0  *  ip 
Sj(b,s)  =  Ij (sb)  [y  +  l0ge  s] 


Similar  formulas  are  obtained  for  the  inside  radius. 


K0(sa)  =  Gn(sa)  -  Mn(sa)  -  Sn(a,s)  , 


GQ(sa) 


0 


O' 


O' 


(8) 


(-m . i-fe). 


(9) 


MQ(sa)  =  IQ(sa)  logg 


SQ(a,s)  =  IQ(sa)  [y  +  loge  s]  , 


and 


^(sa)  =  -  GjCsa)  +  Mj (sa)  +  S^a.s)  , 


(10) 

(11) 


(12) 

(13) 

(14) 

(15) 


(16) 


where 
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(17) 


G]_(s  ,a) 

(sa) 

S]_(a,s) 


1 

(r-1) !r! 


r=l 

!lCsa)  loge  ^ 
I1(sa)  [y  +  log 


1 


1 

r 


(18) 

(19) 


We  now  eliminate  Bessel  functions  from  the  formulas  for  the  cross 
products.  We  obtain 


L0j0(s)  =  I0(sa)  [G0(sb)-M0(sb)]-I0(sb)[G0(sa)-M0(sa)]  (20) 

L0jiCs)  =  I0(sa) [-G1(sb)+M1(sb)]+I1(sb) [GQ(sa) -MQ(sa) ]  (21) 

L1,0(S)  =  !i(sa)  [G0(sb)-M0Csb)]  +  IoCsb)[-Gi(sb)^i(sb)]  (22) 

Ll^l(s)  =  I1(sa)[-G1(sb)+M1(sb)]-I1(sb)[-G1(sa)+M1(sa)]  .  (23) 

The  logarithmic  singularity  in  s  is  eliminated.  The  only  remaining 
singularity  is  a  simple  pole  at  the  origin  due  to  the  functions  M  (sa) 
and  M.(sb).  1 

The  first  few  terms  of  the  ascending  series  for  the  j (s)  functions 
can  be  obtained  by  multiplication  of  the  appropriate  series!  We  find 


'0,0(s) 

=  -  loge(b/a)  +  [-  (a2+b2)loge(h/a)+(b2-a2)Js2/4  +  ... 

(24) 

0,1 (s) 

=  l/(sb)  +  [2b  loge  (b/a)  -  (b2-a2)/b]s/4  +  ... 

(25) 

l,0(s) 

=  1/ (sa)  +  [-2a  loge  (b/a)  +  (b2-a2)/a] s/4  +  ... 

(26) 
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L1,1(S)  =-(b2-a2)/(2ab)+[ab  loge(b/a)-(b4-a4l/(4ab)]s2/4  +  ...  (27) 

Calculation  of  literal  coefficients  for  higher  powers  of  s  requires 
excessive  algebraic  manipulation.  Numerical  values  are  readily  calculated 
from  recursion  formulas,  using  the  first  terms  in  the  above  expansions 
as  starting  values. 

The  recursion  formulas  for  the  coefficients  are  derived  from  a  set 
of  first-order,  linear  differential  equations  satisfied  by  the  cross 
products. 

i 


L  0,0ts)  “  aLl,0(s)  -  bL0,l(s) 

(28) 

L  0,lts)  =  aLl,l(s)  ■  bL0,0<s) 

-  L01(s)/s 

(29) 

L  l,0(s'  =  aLo,o^s-1  '  bLl,l^s^ 

-  Ll,0(s’/a 

(30) 

L'l,ltS)  -  aL0,l<s>  -  bL1.0(s) 

-  2Lj  ^ (s)/s 

(31) 

assume 

00 

L0,oM  ‘n  a" 

n=0 

(32) 

oo 

L0,lts)  °  yi  un  s" 
n=-l 

(33) 

oo 

Ll,0(s>  =2  Vn  S" 

n=-l 

(34) 

oo 

Ll.l(s)  a" 

(35) 

n=0 
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We  substitute  these  series  into  the  preceding  differential  equations, 
and  the  coefficients  of  each  power  of  s  are  set  equal  to  zero.  The 
following  recursion  formulas  are  obtained: 


(n+2)w  -  a  u  .  +  b 

n  n-1 

Vn-1  "  0 

(36) 

(n+2)vn+i  -  a  tn  +  b 

w  =  0 
n 

(37) 

(n+2)un+i  -  a  wn  +  b 

t  =  0 
n 

(38) 

nt  -  av  +  b  u 

n  n-1  n- 

i  = 0 

(39) 

The  indices  are  shifted  for  convenience  in 

computation. 

-  a  u 

n 

+  b  v 

n 

+  (n+3)w  , 
v  *  n+1 

=  0 

(40) 

a  t 

n 

+  (n+23Vi 

+  b  w 

n 

=  0 

(41) 

b  tn  *  (n*2)Vl 

-aw 

n 

=  0 

(42) 

(n+1)‘n+i  *  b  un 

-  a  u 

n 

=  0 

(43) 

We  find  that  Lq^(s)  and  ^(s)  are  even  functions  of  s,  whereas  Lq  ^(s) 
and  Lj  q(s)  are  odd  functions.  The  initial  values  for  the  recursion 
formulas  are  obtained  from  Eqs.  (24  -  27) 


tQ  »  -  loge (b/ a) 

(44) 

u_i  =  1/b 

(45) 

v_i  =  1/a 

(46) 

wQ  =  -  (b2-a2)/ (2ab) 

(47) 
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and 


tj  =  0 


u0  =  0 


V0  =  0 


w^  =  0 


(48) 

(49) 

(50) 

(51) 


The  coefficients  are  listed  in  Tables  1,  2,  and  3  for  several  values 
of  b  and  a=l .  The  series  converge  exponentially  for  all  values  of  s.  The 
convergence  is  very  rapid  when  the  wall  ratio  b/a  is  close  to  unity, 
facilitating  calculations  for  thin-walled  cylinders.  Eqs.  24  thru  27 
can  be  expressed  as 

L0>0Cs)  -  -  l°geCb/a)  *  R00(s) 

L01Cs)  -  1/Csb)  ♦  R01(s) 

Ll,0(s)  =  1/(sa)  +  Rl,0(s) 
bjjCs)  =  -  (b2-a2)/2ab  .  R{  L(s) 

The  graphs  of  the  R  function  for  b  =  1.1,  1.2,  1.25,  1.3,  1.4,  and  1.5 
are  shown  in  Figures  1  thru  12. 

IV.  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS  FOR  MODIFIED 
BESSEL  FUNCTIONS 

We  obtain  asymptotic  expansions  of  the  cross  products  by  expressing 
the  Bessel  functions  of  Eq.  (2)  in  terms  of  Hankel  asymptotic  expansions. 
We  assume 


Im  (s)  >  0 


(52) 


to  be  consistent  with  formulas  for  calculating  stresses  by  the  theory 
of  residues.  Then6 


G.N.  Watson ,  A  Treatise  on  the  Theory  of.  Bessel  FimctionSj  The 
MacMillan  Co.  New  York ,  1948.  See  Eq.  (2)  Page  203 ,  and  the  discussion 
of  Stokes  phenomena  which  follows. 
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TABLE  1.  COEFFICIENTS  FOR  THE  ASCENDING  SERIES  FOR  THE  CROSS  PRODUCTS  OF 
MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIOS  B  =  1.1  AND  B  =  1.2 
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TABLE  2.  COEFFICIENTS  FOR  THE  ASCENDING  SERIES  FOR  THE  CROSS  PRODUCTS  OF 
MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIOS  B  =  1.25  AND  B  =  1.3 
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TABLE  3.  COEFFICIENTS  FOR  THE  ASCENDING  SERIES  FOR  THE  CROSS  PRODUCTS  OF 
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OJ 

X 

O'  OJ 

if) 

X 

pH 

<t  h- 

o 

o 

o 

o  pH 

rH 

pH 

rvi 

OJ  OJ 

1 

1 

1 

1 

1 

1 

1 

i 

1  I 

LU 

LU 

LU 

LU  LU 

LU 

LU 

LU 

LU  LU 

O 

r- 

ro 

st  X 

X 

PH 

OJ 

LO  OJ 

OJ 

o 

X  X 

LO 

ro 

pH  co 

o 

ro 

00 

X  X 

Is- 

o 

OJ 

O'  h- 

o 

<r 

OJ 

<t  -4 

ro 

r- 

o 

st  LO 

z 

X 

X 

x 

O  X 

x 

X 

ro 

CO  pH 

GO 

O' 

<t 

pH  Is- 

pH 

pH 

in 

x  <t 

rH 

in 

o 

O'  o' 

LO 

o 

0s  o 

O' 

in 

X 

OJ  o 

ro 

LO 

LO 

pH  X 

o 

OJ 

h- 

pH  O' 

ph 

O' 

O 

OJ  X 

O' 

Is- 

OJ 

X  X 

a 

X 

pH  X 

• 

• 

• 

•  • 

• 

• 

• 

•  • 

1 

i 

1 

1 

t 

i 

1 

1 

1 

pH 

ro 

LO 

Is-  o 

ro 

X 

O' 

ro  x 

o 

o 

O 

O  pH 

r-H 

pH 

pH 

OJ  OJ 

1 

i 

1 

1  « 

1 

1 

1 

t  1 

LU 

LU 

LU 

LU  LU 

LU 

LU 

LU 

LU  LU 

O' 

O 

ro 

pH  CO 

<t 

if) 

Is-  LT 

X 

O' 

r-H 

OJ  s* 

X 

Is- 

<r 

<t  ro 

pH 

ro 

o 

r-H  OJ 

CO 

CO 

X 

O  O' 

— * 

X 

OJ 

O' 

OJ  <t 

<t 

Is- 

X 

OJ  CO 

z 

X 

X) 

LO  LO 

LO 

r- 

X 

st  OJ 

w 

ro 

X 

O' 

O'  OJ 

X 

ro 

Is- 

X  if) 

> 

X 

r- 

<t 

pH  CO 

OJ 

o 

X  o 

Is- 

<t 

n 

O  X 

o 

ro 

if) 

X  Is- 

pH! 

ro 

o 

<t 

o 

X 

Is- 

if)  X 

Is- 

O' 

<t 

pH  CVJ 

ro 

OJ 

r-H 

X  <f 

r-H 

ro 

X 

X 

o 

ro 

X 

O' 

OJ  X 

o 

o 

o 

o 

r-H 

pH 

pH 

r-H 

OJ  OJ 

1 

l 

1 

i 

1 

1 

1 

1 

i  i 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU  LU 

O' 

if) 

r- 

X 

X 

pH 

r-H 

ro 

<t  O' 

Is- 

o 

X 

<t 

if) 

o 

O 

<t 

pH  O 

Is- 

LO 

Is- 

X 

X 

OJ 

>fr 

r-H 

OJ  <t 

-p**- 

pH 

O' 

X 

o 

OJ 

ro 

ro 

ro 

X  Is- 

z 

X 

O' 

O' 

X 

in 

in 

X 

LO  pH 

w 

r-H 

pH 

CO 

LO 

OJ 

<t 

r- 

x  o 

i 

OJ 

X 

<t 

O' 

CO 

r- 

OJ 

X 

Is- 

OJ 

r- 

pH 

o 

o 

pH 

in 

Is-  X 

Is- 

pH 

.21 

OJ 

r-H 

%t 

<t 

o 

r-H 

X 

r-H 

O' 

pH 

Is- 

r-H 

OJ  X 

pH  X 

1 

1 

1 

1 

1 

1 

1 

1 

1  1 

o 

OJ 

Is- 

O' 

OJ 

LO 

X 

r-H 

<t 

o 

o 

o 

o 

o 

pH 

pH 

pH 

OJ 

OJ 

♦ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

O' 

X 

OJ 

X 

OJ 

X 

X 

LO 

X 

if) 

LO 

X 

OJ 

O' 

X 

r-H 

X 

LO 

in 

in 

Is- 

OJ 

if) 

o 

X 

rvj 

pH 

o 

r— « 

o 

Is- 

O' 

<t 

X 

in 

OJ 

in 

Is- 

ro 

S) 

Is- 

<t 

X 

o 

LO 

<t 

X 

O' 

if) 

X 

X 

X 

<t 

CO 

O' 

h- 

OJ 

O' 

o 

OJ 

ro 

X 

pH 

O' 

OJ 

OJ 

pH 

<t 

ro 

o 

CO 

pH 

in 

o 

OJ 

X 

pH 

OJ 

OJ 

X 

o 

O' 

pH 

pH 

OJ 

r-H 

X 

OJ 

if) 

X 

<t 

CO 

i 


pH 

ro 

in 

Is- 

O 

ro 

X 

O' 

ro 

in 

o 

o 

o 

o 

r-H 

pH 

r-H 

r-H 

Oj 

Oj 

1 

i 

1 

t 

1 

1 

1 

1 

1 

1 

tu 

LU 

LU 

1U 

LU 

LU 

LU 

LU 

UJ 

LU 

pH 

O 

xt 

CO 

X 

r-H 

CO 

pH 

X 

O 

OJ 

in 

o 

o 

CO 

pH 

in 

X 

ro 

pH 

xt 

Is- 

CO 

in 

in 

OJ 

OJ 

X 

o 

pH 

— * 

O' 

0J 

X 

OJ 

xt 

CO 

X 

X 

O' 

X 

z 

O' 

<t 

X 

OJ 

X 

r-H 

O' 

o 

OJ 

pH 

w 

r-H 

in 

X 

CO 

in 

OJ 

o 

in 

<t 

z> 

O 

CO 

o 

O' 

Is- 

xt 

in 

OJ 

r-H 

r- 

r-H 

CO 

o 

xt 

o 

Oj 

Is- 

OJ 

X 

X 

>* 

Is- 

r- 

CO 

O' 

in 

r- 

in 

OJ 

X 

X 

r-H 

ou 

OJ 

OJ 

r-H 

O' 

pH 

r-H 

OJ 

xt 

Is- 

O' 

OJ 

in 

X 

pH 

<t 

o 

o 

o 

o 

o 

r-H 

rH 

r-H 

OJ 

OJ 

I 

I 

1 

( 

1 

1 

1 

1 

1 

1 

LU 

LU 

LU 

UJ 

LU 

LU 

UJ 

UJ 

LU 

LU 

in 

O' 

ro 

o 

X 

O 

in 

pH 

O 

CO 

Is- 

Is- 

if) 

o 

X 

O' 

X 

Is- 

pH 

in 

Is- 

X 

o 

pH 

r-H 

in 

Is- 

0J 

X 

Is- 

0' 

in 

OJ 

in 

xt 

r-H 

>* 

in 

co 

S* 

O' 

O' 

O' 

PH 

r- 

CO 

r-H 

CO 

in 

OJ 

Is- 

xt 

o 

X 

0J 

Oj 

X 

Is- 

X 

o 

pH 

in 

O' 

OJ 

Is- 

Is- 

xf 

OJ 

xt 

X 

o 

<t 

X 

OJ 

in 

in 

o 

Is- 

Is- 

r-H 

pH 

X 

X 

X 

OJ 

O' 

OJ 

r-H 

h- 

OJ 

in 

in 

xt 

CO 

OJ 

<t 

X 

O' 

OJ 

in 

X 

pH 

<t 

X 

o 

o 

o 

o 

fH 

r-H 

r-H 

OJ 

0J 

0J 

1 

1 

1 

i 

1 

« 

1 

1 

1 

1 

UJ 

UJ 

LU 

UJ 

LU 

LU 

LU 

UJ 

LU 

LU 

o 

O' 

Is- 

o 

in 

in 

X 

X 

o 

in 

X 

o 

ro 

X 

in 

X 

Is- 

<t 

r-H 

X 

CO 

OJ 

pH 

O' 

0J 

r- 

X 

X 

• 

in 

X 

X 

OJ 

r-H 

o 

r-H 

X 

ro 

X 

pH 

z 

in 

X 

X 

O' 

o 

O 

pH 

w 

CO 

o 

<t 

X 

Is- 

O 

pH 

O' 

Is- 

II 

h- 

O' 

Is- 

CO 

Is- 

o 

0' 

X 

in 

r-H 

o 

X 

X 

O' 

O' 

Is- 

O' 

CO 

0j 

X 

ro 

O' 

pH 

Is- 

o 

X 

o 

O' 

pH 

X 

pH 

X 

Is- 

0J 

X 

X 

O' 

X 

0J 

pH 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

r—i 

II 

<r 

z 

r-H 

OJ 

ro 

in 

X 

Is- 

X 

O' 

o 

r-H 

CO 

in 

X 

r-H 

CO 

X 

O' 

0J  X 

o 

o 

o 

o 

r-H 

pH 

pH 

r-H 

0J  0J 

1 

1 

1 

1 

1 

1 

1 

1 

1  1 

LU 

UJ 

UJ 

UJ 

LU 

LU 

LU 

LU 

LU  LU 

in 

Is- 

O' 

X 

X 

Is- 

X 

O'  O' 

in 

CO 

OJ 

O' 

X 

Is- 

<t 

in 

0J 

O  X 

o 

X 

in 

X 

0' 

<t 

0J 

ro 

ro  x 

• 

o 

O' 

in 

X 

X 

in 

pH 

in  r- 

pH 

z 

0J 

<t 

in 

X 

O' 

Is- 

0s  0J 

II 

w 

o 

o 

Is- 

ro 

X 

CO 

pH 

X 

o  o 

i— 

o 

0J 

X 

r-H 

X 

CO 

O' 

X  p-» 

CD 

O' 

0J 

X 

X 

Is- 

O' 

o 

<t 

x  r- 

X 

r-H 

0J 

CO 

0' 

in 

O' 

Is- 

0J  X 

pH 

0J 

r-H 

O' 

r-H 

pH 

r-H 

pH  00 

o 

• 

1 

1 

1 

1 

1 

1 

1 

1 

1  1 

pH 

II 

<r 

z 

pH 

2 

3 

4 

5 

6 
7 

X  O'  G 

19 


4017409E-26 


ROO 


It 

jl 


1 


20 


Figure  1.  Graph  of  Rn  n(sl  and  R0  -.(s)  for  Wall  Ratio  =  1.1 
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Figure  2.  Graph  of  R 
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Figure  3.  Graph  of  Rn  nCs)  and  Rn  ,(s)  for  Wall  Ratio 
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Figure  4.  Graph  of R 
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Figure  5.  Graph  of  Rn  n(.$l  and  Rn  ,ts)  for  Wall  Ratio 
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Figure  6.  Graph  of  R nC$l  and  R,  ,(sl  for  Wall  Ratio  =  1.25 
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Figure  7.  Graph  of  Rn  n(_s)  and  Rn  ,(sl  for  Wall  Ratio 
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Figure  8.  Graph  of  R 
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Figure  9.  Graph  of  R 
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Figure  10.  Graph  of  R,  n(jsl  and  R,  ,  (si  for  Wall  Ratio  =1.4 
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Figure  12.  Graph  of  R,  n(_s)  and  R,  .(si  for  Wall  Ratio 


I0(sa)  = 

[eSaPQ(sa)  + 

ie"saq0(sa)] 

//  2irsa 

> 

(53) 

K0(sa)  = 

ire_SaQ0(sa)/ 

liTsa  , 

(54) 

Ij(sa)  = 

[eSaP1(sa)  - 

i  e"saQ1(sa) 

|  //  2irsa  , 

(55) 

Kj(sa)  = 

■sa  A  ^ 

n-e  Qjfsa) 

/  t/  27Tsa  , 

(56) 

where 


PQ(sa)  = 


Q0(sa) 


=  1  + 


V 


P1(sa) 


Q1(sa) 


1!  (8sa) 
2 


1  - 


1 


=  1  - 


=  1  + 


1 ! (8sa) 
1*3 

1 ! (8sa) 


1*  3 


1 ! (8sa) 


2  2 
1  -3Z 

2 ! (8sa) 2 

2  2 

2 ! (8sa) i 
12«  3«  5 

2! (8sa) 1 

l2,  3*  5 

2! (8sa)Z 


(57) 

(58) 

(59) 

(60) 


Similar  formulas  for  PQ(sb),  QQ(sb) ,  P^sb),  and  Q:  (sb)  are  obtained  by 
substituting  b  for  a  in  the  preceding  equations. 


The  following  formulas  for  the  L.  . (s)  functions  are  obtained  by 
successive  substitutions.  1 > ^ 

'P0(sa)Q0(sb)-es(b_a)P0(sb)Q0(sa)J  /  ^4abs2 


J0,0(s)  =  [e"S(b"a)p  f-M..sO>-«), 


■  [■ 

[• 


L1  o(s)  =  |e"s(b‘a)P1(sa)Q0(sb)  +  es(b_a)P()(sb)Q1  (sa)J  /^abs' 


J1 }  1  (s)  =  [e"S<"b"a')p1(sa)Q1 


(sb)  -  es(b"a)P 


l  (sb)qi  (sa)j  /^4abs" 


(61) 


e-s(b"a)Po(sa)Q1(sb)  +  es(b"a)P1(sb)Q0(sa)J/  ^4abs2  (62) 


(63) 

(64) 
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Next,  we  eliminate  the  exponential  functions  from  the  subsequent 
analysis.  Assume  n 


(65) 


then  Uifj(s)  and  V^fs)  will  also  satisfy  Eqs.  (2).  Finally  we  write 


°i.J«  =e"Stb'a)  . 


(66) 


and 


Vi_.Cs)  =  es (k-a)  Fi  j(s) 


(67) 


On  substituting  these  expressions  for  ^.(s)  and  V.  .  (s)  we  obtain  the 

following  sets  of  differential  equations  for  E.  . (s)  and  F.  . fs) . 

i ,  i  i .  i  J  ’ 


F 

0,0 

=  a  El,0 

■  bE0,l  * 

(b-a>  Eo,o  • 

(68) 

EM 

■  a  El,l 

■  bEo,o  * 

[(b-a)- ?]  E0,l  • 

(69) 

"Vo 

‘  a  E0,0 

-  «1,1  * 

[Cb-a)-  ?]  El,0  " 

(70) 

eM 

■  a  Eo,i 

-  bEl,0  * 

[(b-a)- 1]  El,l  • 

(71) 

! 

F 

0,0 

■  a  "l.O 

-  bF0,i  - 

(b-«  F0,0  • 

(72) 

F'o,l 

-  bFo,o  - 

[(»-=>)-  |]f0_i  , 

(73) 
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F  1,0  “  a 


F0,0  "  bFl,l  "  s  ^  Fl,0 


F'i,i  =aFo,i  -bFi,o-  -h  F1,1 


(741 

(75) 


We  assume* 


'0,0 


*n  s-n 


n=l 


J0,l 


u  s 
n 


-n 


n=l 


,10  .  s-»  , 


n=l 


'1,1 


£  "n 


n=l 


We  can  prove  that 


0,0 


£  c-i>n*V" 


n=l 


E,  ..n+l  -n 
(-1)  uns 


n=l 


(76) 


(77) 


(78) 


(79) 


(80) 


(81) 


xo  -  £  c-»n*Vn 


n=l 


(821 


*The  coefficients  t.3  u. ,  v.3  w.  are  new  and  are  not  the  same  as  in 
Eqs.  (32) -(35) .  ^  ^  ^  ^ 
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i,i  ■  E 


(83) 


n=l 


so  a  separate  derivation  for  the  F  .  (s)  functions  is  not  required.  The 
initial  terms  are 


tj  =  l//4ab 

(84) 

Uj  =  l/Aab  , 

(85) 

Vj  =  l//4ab  , 

(86) 

and 


=  l//4ab 


(87) 


The  series  expansions  for  E^,  EQ  ^  E^,  and  Ej  x  are 

substituted  in  Eqs.  (68)  -  (71)  and  the  coefficients  of  each  power  of 

s  set  equal  to  zero.  We  examine  the  coefficients  of  1/s,  1/s2  and  1/s3 
in  detail.  We  find  from  terms  in  1/s 


av1  -  buj  +  (b-a)tx  =  0 


aWj  -  btx  +  (b-a)u  =  0  , 


at1  -  bWj  +  (b-a)v1  =  0  , 

au j  -  bVj^  +  (b-a)w1  =  0 


(88) 

(89) 

(90) 

(91) 
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We  see  these  equations  are  consistent  with  the  initial  conditions  given 
in  Eqs .  (84)  -  (87). 

2 

The  coefficients  of  1/s  lead  to  the  following  equations: 


av2  -  bu2  +  (b-a)t2  =  -  ^  ,  (92) 
aw2  -  bt2  +  (b-a)u2  =  0  ,  (93) 
at2  -  bw2  +  (b-a)v2  =  0  ,  (94) 
au2  -  bv2  +  (b-a)w2  =  w.^  .  (95) 


We  find  these  equations  are  not  linearly  independent,  and  must  be 
supplemented  by  equations  obtained  from  the  coefficients  of  1/s^. 


(b-a)t3 

-bu3 

+av3 

=  -2t2 

(96) 

"bt3 

+  (b-a)u3 

+aw3  =-u2 

(97) 

at3 

+(b-a)v3 

-bw3  =-v2 

(98) 

au3 

-bv3 

+  (b-a)w3  =  0 

(99) 

Consequently, 

2t2  +  U2  +  V2 

=  o  . 

(100) 

We  now  solve  Eqs.  (92),  (93),  (94),  and  (100),  obtaining 

t2  =  t1(b-a)/8ab  ,  (101) 
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u2  =  t1(3a+b)/8ab  f  (102) 

v2~~  t1(a+3b)/8ab  ,  (103) 

W2  =  (3a-3b)/8ab  .  (104) 

We  proceed  with  the  general  ca^e  in  a  similar  manner,  drawing  on 
coefficients  from  1/s  and  1/s  for  a  set  of  linearly  independent 
equations. 

— (n-l)tn_1  (105) 

+  a«n  —  (n-2)un-]_  (106) 

-  b  wn  =-(n-2)vn_1  (107) 

+  (b-a)wn=-(n-3)wn_1(108) 

(-n"1')tn-l  +  (n-2)un_i  +  (n"2)vn_i  +  fn_3^wn_i  =  0  •  (109) 

Now  replace  n  by  (n+1) .  We  find 

n  tn  +  (n-l)un  +  (n-l)vn  +  (n-2)wn  =  0  .  (no) 

To  simplify,  add  Eq.  (105)  to  Eq.  (106),  (107)  and  (108)  in  turn.  We 
find 


(b-a)t  -  bu 
n  n 


+  av 


n 


-  b  t  +  (b-a)u 
n  ^  J  n 


a  t 


n 


+  (b-a)v 


a  u 


n 


-  b  v 


On  adding,  we  find 


-  at 


n 


-  au  +  av  +  aw 
n  n  n 


n-1 


(111) 
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(112) 


bt  -  bu  +  bv  -  bw  =  y  „ 
n  n  n  n  7n-l  9 

(b-a)tn  -  (b-a)un  -  (b-a)vn  +  (b-a)wn  =  zn_1  ,  (113) 

where 

Xn-1  =  '  ^'^Vl  '  (n"2)un-l  »  C114) 

Vl  =  -  (n-1)1:n-l  -  (n-2)vn-l  ’  (US) 

Zn-1  =  '  (n-1)1:n-l  -  (n-3)wn-l  *  (116) 

We  now  solve  Eqs.  (110)  -  (113)  algebraically  by  Cramer* s  rule 
Atn  =  [ (4n-6)a2-  (12n-16)ab  +  (4n-6)b2] [n-l]t  x 

+  [(4n-6)(-ab+b2)][n-2]un_1 

(117) 

+  [(4n-6) (a2-ab)] [n-2]vn_1  +  [ (-4n-4) (ab) ] [n-3]wn  , 

Aun  =  [-(4n-6)(ab-b2)][n-l]tn_1 

+  [ (4n-2)a2  -  (12n-12)ab  +  (4n-6)b2] [n-2]un  , 

(118) 

+  [-(4n-4)ab] [n-2]vn_1 
+  [(4n-2)(a2-ab)][n-3]wn_1  , 
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Avn  =  [(4n-6) (a2-ab)] [n-l]tn_1  +  [- (4n-4)ab] [n-2]un 
+  [(4n-6)a2-(12n-12)ab+(4n-2)b2] [n-2]vn 
+  [-(4n-2) (ab-b2)] [n-3]wn  1  , 

Awn  =  [-(4n-4)ab][n-l]tn_1  +  [ (4n-2) (a2-ab) ] [n-2]u 
+  [- (4n-2) (a2-ab)] [n-2]vn  1 
+[(4n-2)a2-(12n-8)ab+(4n-2)b2] [n-3]wn  1 

where 


(119) 


(120) 


A  =  16(n-l)ab  (b-a) ,  n>l. 


(121) 


These  formulas  were  used  in  an  algorithm  for  the  j  (s)  functions 

which  has  been  incorporated  in  the  thick-walled  cylinder  problem.  The 
variable  s  was  taken  sufficiently  large,  so  the  asymptotic  series  could 
be  truncated  when  a  given  term  was  sufficiently  small,  well  before  the 
smallest  term  of  the  asymptotic  series.  Although  the  formulas  are 
cumbersome,  computations  were  extremely  rapid. 

The  coefficients  are  listed  for  various  wall  ratios  in  Tables  4  thru 
9  and  graphs  of  the  ^ (s)  given  in  Eqs .  61  thru  64  for  positive  real 

s  are  shown  in  Figures  13  thru  24. 

V.  ASYMPTOTIC  SERIES  FOR  SQUARES  AND  PRODUCTS  OF  MODIFIED 
BESSEL  FUNCTIONS 

The  asymptotic  series  for  the  _.  (s)  functions  eliminated  exponential 

overrun  for  all  wall  ratios  of  interest,  but  cancellation  of  leading 
terms  still  occurred,  and  led  to  serious  computational  difficulties  for 
thin-walled  cylinders  when  the  arguments  for  the  Bessel  functions  were 
large.  New  dependent  variables  were  selected  in  which  the  differences 
of  nearly  equal  numbers  were  obtained  by  analysis,  thus  bypassing  the 
numerical  problem.  These  new  variables  were  suggested  by  the  Laplace 
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11153463335E+35  . 10484203123E+34 


TABLE  4.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 
OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.10  (CON 
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TABLE  5.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 
OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.20 
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TABLE  5.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 

OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.20  (CONT'D) 
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TABLE  6.  COEFFICIENTS  FOB  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 
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TABLE  6.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 

OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.25  (CONT'D) 
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TABLE  7.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 
OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.30 
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TABLE  7.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 

OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.30  (CONT'D) 
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TABLE  8.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 

OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.40  (CONT'D) 
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TABLE  9.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 
OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.50 
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TABLE  9.  COEFFICIENTS  FOR  ASYMPTOTIC  EXPANSIONS  OF  CROSS  PRODUCTS 

OF  MODIFIED  BESSEL  FUNCTIONS  FOR  WALL  RATIO  B=  1.50  (CONT»D) 
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Figure  13.  Graph  of  Ln  nCsl  and  LQ  ,(_s)  for  Wall  Ratio 
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Figure  18.  Graph  of  L,  As)  and  L,  ,(s)  for  Wall  Ratio 


LOO  x  LOl  B  =  i . 30 


■ 


o 

* 


o 


CO 


CO 

II 


63 


Figure  19.  Graph  of  Ln  n(_s)  and  L0  -.(s)  for  Wall  Ratio 
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Figure  23.  Graph  of  Ln  nUl  and  Ln  .(si  for  Wall  Ratio 
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expansion  of  the  determinant  in  Eq.  (.1);  we  note  for  instance 


Vp> 

^IjtP) 

ij(p) 

P  I0(P) 

=  P[I02cP)-Ii2(p)]"C2_2v)Il2(P)/P 


Cancellation  occurs  in  the  bracketed  term  when  |p|  is  large. 
Let 

Fj  =  *  [z  Iq2(z)  -  z  I12(z)]  , 

f2  =  tz  I0(*)  K0(z)  +  zIjCz)  Kj (z)  J  , 

F3  !  [zK02(z)  -  zK12Cz)]/»  , 

F4  *•  [Zl0(z)  Kj (z)  ♦  zIjCz)  Kq(z)]  . 

The  last  equation  and  the  Wroqskian  relation 
I0(z)K1(z)  +  I1(z)K0(z)  =  1/z 

give 


(122) 


(123) 

(124) 

(125) 

(126) 


(127) 


F  =  1 

4  A*  (128) 

Derivatives  of  F^,  F2,  and  F^  are  given  by 


=  »  [I02(z)  +  Ij2(z)]  , 

(129) 

=  7T  [4  I0(z)  IjCz)  -  2I12(z)/z]  , 

(130) 

7  — r— — 

Thomas  Muir ,  "A  Treatise  on  the  Theory  of  Determinants".  Dover  Publica¬ 
tions  NY,  1960.  See  p.  94  for  expansion  of  a  4th  order  determinant. 
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[I0(z)  K0(z)  -  IjCz)  Kj (z) ]  , 

(131) 

2  [I1(z)K0(z)-K1(z)I0(z)  +  I1(z)K1(z)/z]  , 

(132) 

[Kq (z)  +  Kj (z)]/ir 

(133) 

-  [4  K0(z)K1(z)+2Kj(z)/z]/tt  . 

(134) 

The  squares  and  products  of  the  individual  Bessel  functions  are  readily 
obtained  from  these  formulas.  We  note  that  cancellation  will  occur  in 
Eqs.  (123),  (125),  (131),  and  (132)  if  |z|  is  large. 

Let 


f=a1f1+a2f2+a3f3+a4f4 


(135) 


then  F  satisfies  the  third  order  differential  equation 

7  IM  2  "  ?  » 

z  F  +  2z  h  -  (4z  +z) F  +  F  =  A4  . 


(136) 


The  functions  F^,  F2,  and  F3  satisfy  the  homogeneous  equation  obtained 

by  setting  A4  =  0;  F.  is  a  particular  solution  of  the  non-homogeneous 
equation. 

Next,  set  sa  =  z  in  Eqs.  (53)  -  (60),  and  define  the  functions 


Vz)  =  po(z)  ’  pi(z) 

(137) 

G2(Z)  =  P0(Z)  Q0(Z)  +  P1(Z)  Ql(z) 

(138) 

G3(z)  =  Qq (z)  -  Qj(z)  , 

(139) 

g4(z)  =  p0(z)Qi(z)  +  P1(z)Qq(z)  . 

(140) 
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It  is  evident  that 


G4(z)  =  1 

on  account  of  Eq.  (127). 

We  also  define  the  functions 

Hj (z)  =  e2zG1(z)  , 
h2Cz)  =  G2(z) 

H3(z)  =  e'2zG3(z)  . 

Then 


(141) 


(142) 

(143) 

(144) 


Fl(z)  =  j  [Hl(z)  ♦  2i  H2(z)  -  H3(z)] 
f2Cz)  =  j  [H2(z)  ♦  iH3(z)]  , 

F3(z)=Ih3(z)  . 


(145) 

(146) 

(147) 


equation6ar  th&t  thC  funCtions  H1 ’  H2’  and  H3  satisfy  the  differential 


3  M  *  2  "  2 

z  H  +  2z  H  -  (4z  +z)  H'  +  H  =  0 


(148) 


since  they  are  linear  combinations  of  F.  ,  F  ,  and  F 

12  3 

Assume 


H 


az 


(149) 
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7  Ml  7  0  M 

z  G  +  (3az  +2z  )G 


+ 


(150.) 


[ (3a2-4)  z3  +  4  a  z2-z]  G* 


+  [(a3-4a)z3  +  (2a2z2)  -  az  +  1]  G  =  0 


Wc  let  a  2,0,  and  -2  in  turn,  obtaining  the  following  differential 
equations : 


3  mi  3  2"  3  2  »  2 

z  Gj  +  (6  z  +2z  )G1  +  (8 z  +  8z  -z)G1  +  (8z  -2z+l)G1 


3  2  "  3 

z  G2  +  2z  G2  +  (4z  +z)G2  +  G2  =  0  , 


0  ,  (151) 

(152) 


7  It  1 

ZG3 

+  (-6z  +12z2)G^+ (-8z3-8z2-z)G^+ (8z2+2z+1)G^  =  0 

(153) 

We  assume 

Gx(z) 

=  al,0  +  al,l  Z_1  +  al,2Z"2  +  ••• 

(154) 

G2(z) 

=  a2,0  +  a2,l  z  1  +  a2,2Z  2  +  • • • 

(155) 

g3(z) 

=  a3,0  +  a3,2  Z  1  +  a3,2Z  2  +  • • • 

(156) 

The  coefficients  are  obtained  by  substituting  each  series  in  the  appropri¬ 
ate  equation,  carrying  out  the  required  algebraic  manipulation,  and  then 
setting  the  coefficients  of  each  power  of  z  equal  to  zero.  The  initial 
terms  are  obtained  from  the  asymptotic  expansions  for  P_(z),  0  (zl  P  fzl 
and  Qj (z) .  We  find  0  Q  J’ 


(8k-8)al ,k" (6k2-14k+6)ai ,k_i+ (k3-5k2+7k-3)ai  k_2  *  0  ,  (157) 
ai,0  =  °»  al,l  =  J  ’  al , 2  "  J  ’  (158) 
4ka2,k-(k3-5k2+7k-3)a2,k_2  =  0  ,  (159) 
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a 


2,0 


(160) 


The  coefficients  of  z  ,  z"3,  z"5  etc.  in  Eq.  (155)  are  zero. 


We  note  that 


G3(z)  =  Gjf-z) 


(161) 


so  that 


(162) 


The  derivatives  of  G^z),  G2(z),  and  G  (z)  are  obtained  by 

differentiating  the  asymptotic  series  in  Eq.  (152),  (153)  and  (154)  term 
by  term.  Derivatives  of  H^z),  H2(z),  and  Hj(z)  are  then  obtained  by 

differentiating  Eqs.  (141)  and  (143).  The  detailed  derivatives  and 
resulting  formulas  will  not  be  given  here. 

Coefficients  for  these  asymptotic  series  are  listed  in  Table  10. 

The  early  coefficients,  which  are  the  most  important,  do  not  increase 
rapidly;  consequently  these  series  offer  an  effective  method  of  computa¬ 
tion.  Characteristic  roots  for  a  thin-walled  cylinder  are  listed  in 
Tables  11  thru  15  for  several  wall  ratios.  These  characteristic  roots 
could  not  be  obtained  either  with  the  original  Bessel  function  subroutine 
or  with  the  asymptotic  series  for  the  L.  . (s)  functions. 

1  y  3 

When  | z |  is  large,  it  is  unnecessary  to  calculate  Fj (z) ,  F_(z),  and 

F3(z),  as  the  characteristic  function  A(s)  can  be  expressed  directly  in 

terms  of  the  P  and  Q  functions  with  exponential  factors.  Combinations 
involving  differences  of  nearly  equal  numbers  can  then  be  expressed  in 
terms  of  the  G  functions.  All  similar  determinants  occurring  in  the 
stress  analysis  of  thick-walled  cylinders  by  Fourier  methods  can  be 
treated  in  a  similar  manner.  First,  replace  the  modified  Bessel  functions 
in  Eq.  (1)  with  the  appropriate  Hankel  asymptotic  series.  These  are  given 
m  Eqs.  (55)  -  (60)  for  the  argument (sa) ;  similar  formulas  hold  for  the 
argument  (sb) .  Next,  eliminate  recessive  terms  like  ie  saQp (sa)//2?rsa 

by  manipulating  columns  of  the  determinant.  We  finally  obtain 
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TABLE  10.  COEFFICIENTS  FOR  THE;  ASYMPTOTIC  SERIES  FOR  THE  SGUAKtS 
AND  PRODUCTS  OF  MODIFIED  HESSll  FUNCTIONS 
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TABLE  10.  COEFFICIENTS  FOR  Thfc  ASYMPTOTIC  SERIES  FO*  THE  SOUARtS 
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TABLE  11.  CHARACTERISTIC  ROOTS  FOR  THIN-WALLED  CYLINDERS,  WALL  RATIO  =  1.010 
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TABLE  12.  CHARACTERISTIC  ROOTS  FOR  THIN-WALLED  CYLINDERS,  WALL  RATIO  =  1.020 
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TABLE  13.  CHARACTERISTIC  ROOTS  FOR  THIN-WALLED  CYLINDERS,  WALL  RATIO  =  1.100 
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TABLE  15.  CHARACTERISTIC  ROOTS  FOR  THIN-WALLED  CYLINDERS,  WALL  RATIO  =  1.250 
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TABLE  15.  CHARACTERISTIC  ROOTS  FOR  THIN-WALLED  CYLINDERS,  WALL  RATIO  =  1.250 
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A(s)  = 


(163) 
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Reduction  of  this  determinant  by  Laplace's  method  leads  to  the  required 
result^ . 


VI.  SUMMARY  AND  CONCLUSIONS 


In  essence,  these  algorithms  permit  precise  calculations  in  an 
extended  region  of  the  complex  s  plane,  so  the  Fourier  integrals  can  be 
evaluated  either  by  quadratures  or  by  residue  theory.  The  ascending 
series  for  the  L^  . (s)  functions  express  A(s)  and  related  determinants 


in  a  form  which  is  free  of  logarithmic  singularities.  Consequently  these 
formulas  are  valid  in  the  entire  finite  part  of  the  complex  s  plane;  the 
cut  along  the  negative  real  axis,  needed  for  the  logarithmic  function,  is 
not  needed.  The  asymptotic  expansions  of  the  L^  . (s)  functions  eliminated 

exponential  overrun  and  also  the  Stokes  phenomena. ^  Finally  the  F^(z) 
functions  eliminated  cancellation  of  leading  terms  which  occurred  for 
large  | s | . 


It  is  clear  that  accurate  values  of  the  Bessel  functions  do  not 
ensure  accuracy  in  evaluating  the  Fourier  integrals  arising  in  thick- 
walled  cylinder  theory  when  js|  is  very  large  or  very  small;  special 
attention  for  various  combinations  of  Bessel  functions  is  also  required. 
The  obvious  alternative,  the  use  of  multiple  precision  calculations,  is 
not  standard  practice  at  BRL  and  consequently  was  not  considered. 


O 

A.S.  Elder,  J.N.  Walbert ,  K.  Zimmerman,  "Stresses  near  a  Diseontinuity 
of  Loading  in  Thick  and  Thin  Walled  Cylinders" ,  BRL  Report  in 
preparation. 
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